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Abstracts: A simple method for an unsteady aerodynamic model tuning is proposed in this 
study. This method is based on the direct modification of the aerodynamic influence 
coefficient matrices. The aerostructures test wing 2 flight-test data is used to demonstrate the 
proposed model tuning method. The flutter speed margin computed using only the test 
validated structural dynamic model can be improved using the additional unsteady 
aerodynamic model tuning, and then the flutter speed margin requirement of 15 % in military 
specifications can apply towards the test validated aeroelastic model. In this study, unsteady 
aerodynamic model tunings are performed at two time invariant flight conditions, at Mach 
numbers of 0.390 and 0.456. When the Mach number for the unsteady model tuning 
approaches to the measured fluttering Mach number, 0.502, at the flight altitude of 9,837 ft, 
the estimated flutter speed is approached to the measured flutter speed at this altitude. The 
minimum flutter speed difference between the estimated and measured flutter speed is -.14 %. 

1 INTRODUCTION 

The primary objective of this study is to reduce uncertainties in the unsteady aerodynamic 
model of an aircraft to increase the safety of flight. To this end, a new flutter analysis 
procedure using the validated aeroelastic model is proposed, and the block diagram of this 
new procedure is shown in figure 1 . 



(a) Previous procedure 
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(b) Proposed new procedure 

Figure 1 : Flutter analysis procedure at NASA Dry den Flight Research Center 


Significant efforts [1, 2, 3] have been made in developing corrections to linear aerodynamic 
models to improve correlation with steady-state wind tunnel and flight test data. There has 
been a limited amount of effort in the correction of unsteady aerodynamics for 
aeroservoelastic applications, which has been relatively sparse and ad-hoc when compared to 
the steady-state work that has been performed. 



Figure 2: Aerostructures Test Wing 2 mounted on F15B pylon for the flight testing. 
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The Aerostructures Test Wing (ATW) 2 test article, shown in figure 2, was developed and 
flown at the NASA Dry den Flight Research Center (DFRC) on the F15B test bed aircraft on 
December 15, 2009. To support the envelope expansion, a Test- Validated Finite Element 
Model (TVFEM) was used for the flutter analysis of the ATW2. Flutter boundaries of the 
ATW2, before and after the structural dynamic model tuning [4], are compared with the flight 
envelopes as shown in figure 3 [5]. In this figure, the solid line bounds the ATW2 test 
envelope that is planned for flight, and the dashed line is the 15% margin of the ATW2 test 
envelope. This 15% margin line was designed to match the numerical flutter boundaries 
computed using the TVFEM with the 3% structural damping, the solid line with circular 
marker. The solid line with the diamond marker represents flutter boundaries using the 
TVFEM with the measured structural damping [5], The measured flutter point of the ATW2 is 
also shown in figure 3, using the x marker. 


Flutter Boundaries 

Flutter Margins 

Measured/ 1 . 1 5 (= Vd) 

0% 

Measured (= 1.15 Vd) 

15 % 

TVFEM; after model tuning with measured damping 

32% 

TVFEM; after model tuning with 3% structural damping 

41 % 


Table 1 : Required flutter margins for different types of structural dynamic models 


■: Test Validated FE Model (3% structural damping) 
Test Validated FE Model (measured damping) 

X : Measured Flutter Point 
+ : Measured Flutter Point/1. 15 
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Figure 3: Flutter boundaries before and after structural dynamic model tuning. 


3 


Required flutter margins for the safety of flight were computed and summarized in table 1 and 
figure 3. It should be noted that the series of Ground Vibration Tests (GVTs) and structural 
dynamic model tuning has been performed [5, 6, 7] resulting in the computed flutter 
boundaries, based on the ATW2 configuration in figure 2 and the corresponding TVFEM, 
shown in figure 3. The GVT for the final ATW2 configuration was performed while the 
ATW2 was mounted to the Flight Test Fixture (FTF) [8] in the FTF ground handling cart as 
shown in figure 4. The FTF was sufficiently massive when compared to the ATW2 so that 
cantilevered boundary conditions were used. It may conclude from table 1 and figure 3 that 
when only the structural dynamic model is validated with respect to GVT data, the flutter 
margin required for the flutter certification of the ATW2 should be approximately 40%. In 
addition to the historically stand-alone structural dynamic model, the unsteady aerodynamic 
model should also be validated with the test data to use the 15% flutter margin in the military 
specification [9], 



Figure 4: Flight test fixture and ATW2 on ground handling cart. 


2 UNSTEADY AERODYNAMIC MODEL TUNING PROCEDURE 

A simple technique has been proposed and developed to update unsteady aerodynamic 
models. The technique is based on matching the measured and numerical aeroelastic 
frequencies of an aircraft structure. In defining the optimization problem to match the 
measured aeroelastic frequencies, the variation of the unsteady aerodynamic force was 
selected as the design parameter. The unsteady aerodynamic force is a function of Mach 
number, reduced frequency, and dynamic pressure; which can be obtained based on any 
aerodynamic model. ZAERO [10] code is used in this study. If the Mach number is constant, 
the reduced frequency and dynamic pressure become variables for changing the unsteady 
aerodynamic force. 

Supporting the Aeronautics Research Mission Directorate guidelines, NASA DFRC has 
developed an Object-Oriented Optimization (O 3 ) tool [11] which leverages existing tools and 
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practices, and allows the easy integration and adoption of new state-of-the-art software. 
Unsteady aerodynamic model tuning used in this study is based on the minimization of the 
discrepancies in numerical and measured aeroelastic frequencies. A computer code for 
unsteady aerodynamic model tuning has been developed using the O tool together with the 
pre-processor, ZAERO, and post-processor codes shown in figure 5. 



Frequencies 
Measured from 
Flight Test 


Figure 5: Unsteady aerodynamic model tuning using object-oriented optimization 


2.1 Pre-processor code 

This code reads in design variables generated by the O tool, and then reads modal 
Aerodynamic Influence Coefficient (AIC) matrices that were computed and saved using 
ZAERO code. Modified modal AIC matrices are then created as shown in figure 6. 



Figure 6: Flow-chart of the pre-processor 
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Design variables in this unsteady aerodynamic model tuning are scaling factors for each 
element in the AIC matrices. The AIC matrix, A, at a reduced frequency can be written as: 
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Where m is number of degrees of freedom, n is number of modes, and ay and by are the i-th 
row and j-th column element of the real and imaginary part of the matrix A, respectively. 
Design variables, en,en, ...,e 2 i,e 2 2 , f 11 / 12 , ... J 21 J 22 , -fmn, are defined as: 
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Where ey is the design variable for ay and fy is the design variable for by. The following 
design variable linking options are available for the unsteady aerodynamic model tuning. 


Option 1 : single design variable 

d <2 1 / 6 ] 2 • • • £ mn~fll~fl2 •• • ~fmn 

Option 2: two design variables 

di =en=ei 2 =...=e m n, real part 
d 2 =fn=fi 2 =-=fmn ; imaginary part 

Option 3: columnwise the same design variables (total n design variables) 

di =en=e 2 i = ...=e m i= fn=f 2 i= ... =f m i 

d 2 =ei 2 =ej 2 =... =e m 2 = f[2 = f22 =■■■ =fm2 
dn £ln (?2n • • • £mn~fln~f2n ••• ~fmn i 

Option 4: columnwise the same design variables (total 2 n design variables) 

d\=eu=e2i= ...=e m i 

d2=ei2=e22=... =e m 2 

d n £jn £ 2 n ••• £mni real parts 
d n +l =fll =fn = ••• =fml 

dn+2 = fl2 = f22 = • • • = fm2 

d 2n =fin=f2n=-=fmn, imaginary parts 


Option 5: No design variable linking; total 2 mn design variables. 
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2.2 ZAERO flutter analysis 

Flutter analyses in this study are based on ZAERO code. This code acquires the modified 
modal AIC matrices and performs the matched flutter analysis as shown in figure 7. This 
computer simulation requires the natural frequencies and mode shapes of the aircraft, and in 
this study these modal data are computed using MSC/NASTRAN code [12]. The V-g and V-f 
data are computed and saved for the next post-processing step. Here, V is aircraft speed in 
Keas, g is structural damping, and f is aeroelastic frequency in Hz. 



Figure 7: Flow-chart of the flutter analysis procedure 


2.3 Post-processor code 

This program reads in the V-g and V-f data and the target altitude where the flight test was 
performed. Based on the velocity information V in the V-g and V-f data, corresponding 
altitudes at fixed Mach numbers are computed. Numerical aeroelastic frequencies are 
computed from target and computed altitudes using the cubic splining procedure. 

This program also computes the frequency difference between the numerical and measured 
aeroelastic frequencies. Frequency difference will be an objective function, which will be 
minimized through the use of the O 3 tool. The proposed tuning technique is an unconstrained 
optimization problem that can be solved using a gradient based optimizer [13], a genetic 
algorithm [14], or a big-bang-big-crunch algorithm [15, 16, 17]. A flow-chart of this post- 
processor code is shown in figure 8. 



Figure 8: Flow-chart of the post-processor 
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3 APPLICATION 

During the flight test of the ATW2, a classical bending and torsion type of flutter, as shown in 
figure 9, was observed near a Mach number of 0.502 and a flight altitude of 9,837ft. 



Figure 9: Classical bending and torsion flutter during flight test 


Measured Mach number, flight altitude, and acceleration at the middle and leading edge of the 
wing tip boom after take-off, near flutter, and during flutter are shown in figures 10 and 11, 
respectively. Measured aeroelastic frequencies during the flight test as well as natural 
frequencies computed using the TVFEM and the measured during GVT are shown in table 2. 
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83720 83730 83740 83750 



Time (sec) 


Time (sec) 


(a) Mach number 


(b) Flight altitude 




Time (sec) 


Time (sec) 


(c) Acceleration at middle of boom (d) Acceleration at leading edge of boom 

Figure 10: Measured Mach number, flight altitude, and acceleration at the middle and the leading edge 
of the boom after take off (time steps 83720 sec to 83750 sec) 
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The first model tuning is performed using flight data after take-off, from 83,742 s to 83,745 s 
as shown in figure 10. The ATW2 in this time period is a time invariant system. The average 
flight Mach number and altitude during this time period were 0.39 and 9,934 ft, respectively. 
Unfortunately, measured acceleration data in figure 10 was noisy, and it was quite difficult to 
estimate the first and third aeroelastic frequencies because of high aerodynamic damping. The 
second measured aeroelastic frequency is 40.45 Hz as shown in table 2. An initial aeroelastic 
frequency of 41.12 Hz is computed using ZAERO code as shown in table 3. 


Mode Natural Frequencies Measured Aeroelastic Frequencies 


TVFEM GVT After take off* Near flutter** During flutter*** 


1 

17.45 

17.45(0.623)1 




2 

43.48 

43.72(0.610) 

40.45 

38.99 

37.69 

3 

82.98 

83.66(0.778) 




4 

133.6 

N/A 




5 

153.8 

142.3(0.674) 





*: Time steps 83,742 sec to 83,745 sec (Mach=0.39 & altitude=9,934 ft; time invariant) 

**: Time steps 84,668 sec to 84,672 sec (Mach=0.456 & altitude=9,858 ft; time invariant) 

***: Time steps 84,683 sec to 84,684 sec (Mach=0.502 & altitude=9,837 ft; time varying) 

( )f : Measured damping (%) 

Table 2: Numerical and measured frequencies (Hz) of the ATW2 during the flight test 


Mach 

Measured 

Altitude 

Before Tuning 

Scaling Factor 

After Tuning 

Number 

(Hz) 

(ft) 

(Hz) 

(design variable) 

(Hz) 

0.390 

40.45 

9,934 

41.12 

1.2579 

40.45 

0.456 

38.99 

9,858 

40.10 

1.2719 

38.99 


Table 3: The second aeroelastic frequency before and after unsteady aerodynamic model tuning 

and corresponding scaling factors 

In this model tuning procedure, the aeroelastic frequency difference in the second mode is 
minimized using the design variable linking option 1. The number of target aeroelastic 
frequencies to be matched is one, and therefore the simplest option is selected. In other words, 
there is an unconstrained optimization problem with a single design variable. After model 
tuning, the second aeroelastic frequency of 41.12 Hz becomes 40.45 Hz, and the 
corresponding scaling factor (single design variable) is 1.2579 as shown in table 3. 


Scaling Flutter Speed 


Flutter Frequency 


Comment 


Factor 

(design 

variable) 


Keas 


% 

difference 


Altitude 

(ft) 


Hz 


% 

difference 


Measured 

N/A 

276.4 

0.00 

9,836.9 

37.69 

0.00 

Before tuning 

1.0 

311.3 

13.0 

3,561.5 

37.67 

-0.05 

Use M=0.39 Aero 

1.2579 

277.3 

0.33 

9,670.0 

37.69 

0.00 

Use M=0.456 Aero 

1.2719 

276.0 

-0.14 

9,912.5 

37.68 

-0.03 


Table 4: Measured and computed flutter boundaries at Mach = 0.502 

The saved AIC matrices at Mach = 0.502 is updated using the scaling factor of 1.2579, and 
the updated flutter boundary at this Mach number is summarized in table 4 and figure 12. The 
tuned flutter speed, corresponding altitude, and flutter frequency are 277.3 Keas, 9,670 ft, and 
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37.69 Hz, respectively. It should be noted in table 4 that flutter speed and frequency 
difference after the unsteady aerodynamic model tuning are 0.33% and 0.00%, respectively. 


0.60 

0.55 

0.50 

I °- 45 
§ 

0.40 

0.35 


0.30 
84660 































I 


























I 


z 

z 

z 

z 

z 

z 

z 

_ 

z 

z 

z 

z 

z 

Z I 

IZ 





z 

~ 

z 

z 

z 

z 





— 

— 

— 

— 

— 

— 

— 

- 


— 

— 

— 

— 

— 

— 

— 


— 

— 

— 

— 

— 

— 

— 

— 

— 

















































Z 







_ 





_ 















- 


■H 




— 

— 

— 


— 

— 

— 

— 

- 

- 








— 

- 

— 

r 

— 


— 

— 

— 

— 

— 

— 


1 



— 




1 

V 



5 





-I 

During r lutter 

— 

- 

— 

— 


— 

— 

- 

- 

— 

- 

- 

- 

- 

- 

- 

— 




— 

L 


— 

— 

— 

— 

— 

— 

J 























_ 

- 

— 

_ 


_ 

z 



z 

z 

z 


-- 

- 

- 

- 

— 

— 

- 

- 

- 

- 




























































84670 84680 

Time (sec) 


10000 

9900 

9800 

9700 

9600 

9500 

9400 

9300 

9200 

9100 

9000 


84690 


Near Flutter 


-H fe- 


During Flutter 


84660 84670 84680 

Time (sec) 


84690 


(a) Mach number 
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(c) Acceleration at middle of boom 


84690 


•-c 
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Time (sec) 


84690 


(d) Acceleration at leading edge of boom 


Figure 11: Measured Mach number, flight altitude, and acceleration at the middle and the leading edge 
of the boom just before and during flutter (time steps 84660 sec to 84690 sec) 


The second model tuning is performed using another time invariant system, between time 
steps of 84,660 s and 84,690 s. Relatively flat time histories of flight Mach numbers and 
altitude, between time steps of 84,668 s and 84,672 s, are observed as shown in figure 11. In 
this time period, the average flight Mach number and altitude were 0.456 and 9,858 ft, 
respectively. 


The second measured aeroelastic frequency in these time steps is 38.99 Hz as shown in table 
2. The corresponding aeroelastic frequency computed from the ZAERO simulation with Mach 
0.456 aerodynamics is 40.10 Hz as shown in table 3. Unsteady aerodynamic model tuning is 
performed using these two numerical and measured frequencies, and a scaling factor of 
1.2719 based on the design variable linking option 1, which is also given in table 3. 


Flutter analysis at Mach 0.502 is performed using this new scaling factor, and the flutter speed 
and frequency are summarized in table 4 and figure 12. Flutter speed difference of 13% before 
the unsteady aerodynamic model tuning becomes -0.14% after tuning. 
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Figure 12: Flutter boundaries at Mach=0.502 before and after unsteady aerodynamic model tuning. 

In case of the ATW2, computation time required for completing an unsteady aerodynamic 
model tuning based on option 1 was less than 7 min. Once the scaling factor (design variable) 
is computed, an additional 1 or 2 min is needed for the Fast Fourier Transformation, one more 
flutter analysis at a higher Mach number to compute the updated V-g and V- f data, and 
automatic computations of updated flutter speed and frequency. Therefore, less than 9 min are 
enough to predict more accurate flutter speed based on the current flight test data. 

4 CONCLUSION 

A simple unsteady aerodynamic model tuning based on the direct AIC modification is 
proposed in this study. The value of the unsteady aerodynamic model tuning procedure has 
been shown with the application to the ATW2 flight test data. 

Flutter boundaries and the ATW2 flight test envelope were computed using the TVFEM [5, 
6]. The flutter margins required for the safety of flight were approximately 40% when only 
the structural dynamic model was validated. Excellent flutter speed matching is accomplished 
when the simple unsteady aerodynamic model tuning is applied resulting in flutter speed 
differences of 0.33% and -0.14%. The flutter margin requirement of 15% in the military 
specification can now be used with the test validated aeroelastic model, that is test validated 
structural dynamic and unsteady aerodynamic models. The modeling uncertainties associated 
with the unsteady aerodynamics can be easily minimized through the use of the simple model 
tuning procedure proposed in this study. 

Unsteady aerodynamic model tunings are performed at two time invariant flight conditions, at 
Mach numbers of 0.390 and 0.456. When the Mach number for the unsteady aerodynamic 
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model tuning approaches to the measured fluttering Mach number, 0.502, at the flight altitude 
of 9,837 ft, the estimated flutter speed is approached to the measured flutter speed at this 
altitude. Therefore, we may conclude that the Mach number selected for the unsteady 
aerodynamic model tuning is closer to the measured fluttering Mach number at the same flight 
altitude, and we may get a more accurate scaling factor for the precise flutter prediction. 
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Unsteady Aerodynamic Model Tuning 
for Precise Flutter Prediction 


Chan-gi Pok, Ph.D. 
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NASA Dryden Flight Research Center 



Objectives 



□ The primary objective of this study is to reduce uncertainties in the unsteady aerodynamic model 
of an aircraft to increase the safety of flight. 


□ This model tuning technique is applied to improve the flutter prediction of the Aerostructures Test 
Wing 2. 


□ This work is supported by the Aeronautics Research Mission Directorate (ARMD) Subsonic Fixed 
Wing (SFW) and Supersonics (SUP) projects under Fundamental Aeronautics (FA) program. 
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Flutter Analysis Procedure @ NASA Dryden 


□ Everyone believes the test data except for the experimentalist, and no one believes the 
finite element model except for the analyst. 


❖ 


Some of the discrepancies come from analytical Finite Element modeling 
uncertainties, noise in the test results, and/or inadequate sensor and actuator 
locations. Not the same orientation for each sensor. 


Structural Dynamic 
Finite Element Model 


Weight, C.G., Moment of 
inertia, & GVT data 


□ Flutter Analysis 

❖ Uncertainties in the structural dynamic model are minimized through the use of 
"model tuning technique" 

❖ Based on analytical modes 

□ Validate Structural Dynamic Finite Element Model using Test Data and Update if needed 


❖ 


Use MDAO (Multidisciplinary Design, Analysis, and Optimization) tool with Model 
Tuning Capability or Standalone Model Tuning Code 

> Model tuning is based on optimization. 

■S Design Variables 

• Structural sizing information: Thickness, cross sectional area, 
area moment of inertia, etc. 

• Point properties: lumped mass, spring constant, etc. 

• Material properties: density. Young's modulus, etc. 

S Constraints 


Structural Dynamic 
Model Tuning 


Validated Structural 
Dynamic Model 


Create Unsteady 
Aerodynamic Model 
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Pre-Test Modal & Flutter Analyses of ATW2 



Structural Dynamic & Unsteady Aerodynamic Models 



Cross Sectional Shape: 
NACA-65A004 


□ Structural Dynamic Finite Element Model 

♦> Based on MSC/NASTRAN code 

Use ATW1 Structural Dynamic Finite Element Model (265 nodes) 
> ATW1 & ATW2: Based on same drawing 
❖ Use 10 modes for the flutter analysis 


ining Points 


□ Unsteady Aerodynamic Model 

❖ Based on ZAERO code 

❖ 416 elements 



❖ Select 16 reduced frequencies between 0 & 1 

❖ Mach = .60, .75, .82, and .95 

❖ Linear Theory 

00 A Use Matched Flutter Analysis 
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Splined Mode Shapes 



Mode 1: 17.60 Hz 
First Bending 


Mode 2: 23.26 Hz 
First Torsion 


Mode 3: 93.99 Hz 
Second Bending 
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Frequency (Hz) Damping 


V-g and V-f Curves at Mach = 0.82 Before Model Tuning 



Speed (Keas) 



J 

HH 

3 

F0q 

r* 

he: 


































VA 


A J VV 



















' | 

|V C| V | 















v“* "W- — 

» — 

| y, | x 1 



















T ' 


A |a] 



















flj— p=H — 



i- 


A-j 

p b 


r. ' 

r* 




































































































































nil 



5= 

3=f 

—LI 
— [} 
[3 




r 

r 

H 



* 


ftijj 



j§fc 




*— 

♦ — 


























200 300 400 500 600 700 
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Damping 


Mode 

Frequency 

Modal Participation Factor 

1 

17.60 Hz 

75.0% 

96.6% 

2 

23.26 Hz 

16.8% 

3 

93.99 Hz 

4.8% 

4 

135.4 Hz 

0.0% 

3.4% 

5 

163.1 Hz 

2.6% 

6 

174.5 Hz 

0.0% 

7 

257.5 Hz 

0.5% 

8 

391.6 Hz 

0.0% 

9 

394.3 Hz 

0.1% 

10 

445.6 Hz 

0.3% 


Flutter Mode 

Speed 

Frequency 

Altitude 

1 

407.4 Keas 

22.86 Hz 

15010 ft 
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Finite Element Model Tuning 
of the Aerostructures Test Wing 2 
Using Ground Vibration Test Data 



Structural Dynamic Model Tuning using GVT Data 


□ MIL-STD-1540C Section 6.2.10 

♦♦♦ Test Requirements for Launch, Upper-Stage, & Space Vehicles 

♦♦♦ Less than 3% and 10% frequency errors for the pnman and secondary, modes, respectively 
❖ Less than 10% off-diagonal terms in orthonormalized mass matrix 


□ AFFTC-TIH-90-001 (Structures Flight Test Handbook) 

❖ If measured mode shapes are going to be associated with a finite element model of the 
structure, it will probably need, to be adjusted to match the lumped mass modeling of the 

analysis. 

❖ Based on the measured mode shape matrix ® and the analytical mass matrix M, the following 
operation is performed: 

® t m® 

♦♦♦ The results is near diagonalization of the resulting matrix with values close to 1 on the diagonal 
and values close to zero in the off-diagonal terms. Experimental reality dictates that the data will 
not produce exact unity or null values, so 10 percent of these targets are accepted as good, 
orthogonality and the data can be confidently correlated with the finite element model. 
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Structural Dynamic Model Tuning Procedure 


□ Pak, C.-G., "Finite Element Model Tuning using Measured Mass Properties and Ground 
Vibration Test Data/' ASME Journal of Vibration and Acoustics, Vol. 131, Issue 1, 
February 2009. 

♦> Applied to beam finite element model using in-house FEM 


□ Pak, C.-G. and Shun-fat Lung, "Reduced uncertainties in the flutter analysis of the 
aerostructures test wing," NASATM 2011-216421, 2011. 

♦> Using in-house object-oriented optimization tool 

♦♦♦ Structural model tuning is based on MSC/NASTRAN code. 

> Single performance index for frequency differences 

> Single performance index for off-diagonal terms of the orthonormalized mass matrix 


□ Pak, C.-G. and Shun-fat Lung, "Flutter Analysis of the Aerostructures Test Wing with Test 
Validated Structural Dynamic Model," Journal of Aircraft (accepted for publication 25 
Feb. 2011). 

♦> Extended methods are used 

> Multiple performance indices for frequency differences 

> Multiple performance indices for off-diagonal terms of the orthonormalized mass 
matrix 
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Test Setup: #1 GVT (Strong Back Mounting) 



Sensor: Sticker 


□ Before installing Teardrop Accelerometers & TAO Pressure Sensors 

□ Strong Back Mounting @ Flight Loads Lab 

□ Use Photogrammetry Optical Measuring System 


23 d 

t 

Excitation point 
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Additional Sensors for Flight Test 


Bottom View 



Mounting Panel (Not Shown Fully) 


Strain Bridges (10) 


Piezoelectric Excitation Patch Pairs 


Strain Bridges (4) 


Teardrop Accelerometers 
(2: Epoxied onto Surface) 


TAO Pressure Sensors 


Teardrop Accelerometers (2: Inside Boom) 
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Test Setup: #2 GVT (FTF Mounting) 



□ 

□ 

□ 


After installing 4 Teardrop Accelerometers & 2 TAO Pressure Sensors 
Flightiest Fixture Mounting @ F-15B Hanger 
Flight Test Fixture was lot heavier than ATW2. 

❖ FTF * 500 lb vs. ATW2 = 2.66 lb: 500/2.66=188 » 10 
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Test Setup: #3 GVT (FTF Mounting under F-15) 



& ASA 


Center Fuselage Pylon 




□ After installing under the Center Fuselage Pylon 


□ Add Flexibilities between Flight 
Center Fuselage Pylon 


Br Fuselage Pylon 
t Test Fixture & 


Structural Dynamics Group 


Chan-gi Pak-14 








Results (Frequency Comparisons) 


Mode 

GVT (Hz) 

Before Tuning 

After Tuning (Target #2 GVT data) 

MIL- 

STD 

(%) 

#1: 

Strong 

Back 

#2: 

FTF 

#3: 
FTF & 
F15B 

Freq 

(Hz) 

Error (%) 

Freq 

(Hz) 

Error (%) 

Wrt 

#1 

Wrt 

#2 

Wrt 

#3 

Wrt 

#1 

Wrt 

#2 

Wrt 

#3 

1 

17.24 

17.45 

17.42 

17.60 

2.09 

0.86 

1.03 

17.45 

1.22 

0.00 

0.17 

3 

2 

44.10 

43.72 

43.73 

23.26 

- 47.3 

- 46.8 

- 46.8 

43.48 

- 1.41 

- 0.55 

- 0.57 

3 

3 

84.00 

83.66 

84.14 

93.99 

11.9 

12.4 

11.7 

82.98 

-1.21 

-0.81 

-1.38 

3 

4 

N/A 

N/A 

N/A 

135.4 

N/A 

N/A 

N/A 

133.6 

N/A 

N/A 

N/A 


5 

N/A 

142.3 

143.0 

163.1 

N/A 

14.6 

14.1 

153.8 

N/A 

8.08 

7.55 

10 
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Results (Total Weight, Orthogonality, & MAC) 



Measured 

Before Tuning 

After Tuning 

Total 

Weight 

2.66 lb 

1.76 lb (error 34%) 

2.85 lb (error 7.1%) 

Orthonormalized 
Mass Matrix 


1 

/-r viol 

ate JS 

1 

T Sat 

isfy ^3 

1 

1 

\24.9% 

38.0%^ 

1 

\l.92% 

-4.46%\ 

2 

-.249 


-66A%y 

-.0192 


6.16%^/ 

3 

.380 

-.661 

/ 1 

-.0446 

.0616 

/ 1 

MAC 

Mode 1 

.97 / 

.99 / 

Mode 2 

.70 / 

.99 / 

Mode 3 

,75 / 

.98 / 


MIL-STD & AFFTC-TIH-90-001 Requirements: 10% 
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Flutter Analysis Using Validated Structural 
Dynamic Finite Element Model 



Frequency (Hz) Damping 


V-g and V-f Curves at Mach = 0.82 After Model Tuning 


0.1 
0.0 
- 0.1 
,- 0.2 
-0.3 
-0.4 
-0.5 

200 300 400 500 600 700 


Damping 


Mode 

Frequency 

Modal Participation Factor 

1 

17.45 Hz 

10.7% 


2 

43.72 Hz 

87.2% 

99.5% 

3 

83.66 Hz 

1.6% 
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Speed (Keas) 



Flutter Mode 

Speed 

Frequency 

Altitude 

1 

341.5 Keas 

34.59 Hz 

23475 ft 


Speed (Keas) 
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Altitude (ft) 


Flutter Boundaries vs. Flight Envelope 



: Flutter Boundary Before Model Tuning (3% structural damping) 

: Flutter Boundary After Model Tuning (3% structural damping; FEM based on #1 GVT data) 

: Flutter Boundary After Model Tuning (3% structural damping; FEM based on #2 GVT data) 

: Flutter Boundary After Model Tuning (measured structural damping; FEM based on #2 GVT data) 



Mach 
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Flight Test & Summary of Flutter Margins 




ATW2 Flight Test 



Flight Test 

ryden Flight 
rch Center 

er 15 r 2000 


Flight T est 

ryd e n Flight 
rch Center 

er 1 5, 2009 



Altitude (ft) 


50000 


40000 


30000 


20000 


10000 


0 


-10000 


Summary of Flutter Margins 



: Flutter Boundary Before Model Tuning (3% structural damping) 

: Flutter Boundary After Model Tuning (3% structural damping; FEM based on #1 GVT data) 

: Flutter Boundary After Model Tuning (3% structural damping; FEM based on #2 GVT data) 

: Flutter Boundary After Model Tuning (measured structural damping; FEM based on #2 GVT data) 
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Summary of Flutter Margins (continued) 


Flutter Boundaries 

Flutter Speed Differences 

Measured/1.15 = V d 

0% 

Measured = 1.15 V d 

15% 

Test validated FEM; using #2 GVT data; with measured damping 

32% 

Test validated FEM; using #2 GVT data; with 3% structural damping 

41% 

Test validated FEM; using #1 GVT data; with 3% structural damping 

52% 

FEM; before model tuning; with 3% structural damping 

74% 


Validated Structural 
Dynamic Model 

Validated Unsteady 
Aerodynamic Model 

Recommended* 
Flutter Margins 

ATW2 Case 

Yes 

Yes 

15% 


Yes 

No 

49% 

32 - 52 % 

No 

No 

54% 

74% 


*: POP's JSSGz2006 Guidelines for Flutter Speed Clearance; Faustino Zapata, AFPC May 22-23, 2008 
JSSG(Joint Service Specification Guide) 

Puring Pesign Phase (paper plane phase) : 

□ Can we use the 15% flutter margin??? 

□ If not, then what percentage of flutter margin should be used. 54% ?? 
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Unsteady Aerodynamic Model Tuning 



New Flutter Analysis Procedure @ NASA Dryden 
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Unsteady Aerodynamic Model Tuning 


□ Optimization Problem Statement 

♦♦♦ Objective Function: 

Minimize J = measured aeroelastic frequency - computed aeroelastic frequency 
♦> Design Variables: e^ & f Vj 


e \ 1 a l 1 

e n a n 

— e \n a \n 


fn b n 

f\2 b \2 

i 

s 

s 


a li 

a 12 

— a \n 


~b u 

b 12 

1 

R 

e 2l a 2\ 

6^22^22 

— e 2n a 2n 

+ i 

fiAi 

f22 b 22 

••• f2n b 2n 

A = 

a 2l 

a 22 

— a 2^z 

+ i 

b 21 

b 22 

••• b 2n 

1 

aP 

ft 

3 

ml ml 

e mn a mn _ 


_f mAm\ 

f m2 b m2 

*'* -f mn b mn _ 


a ,n\ 

a m 2 

^mn _ 


Ai 

b m 2 

••• b mn _ 


♦♦♦ Design Variable Linking 

> Option 1: single design variable 

d =e il =e i2 = ’” =e mn = fll = fl2 = - = fmn 


> Option 2: two design variables 

d i =e ii =e i2=- =e mn > real P art d 2 =f u=f i 2 =”-=f mn '> imaginary part 

> Option 3: columnwise the same design variables (total n design variables) 

d l =e ll =e 21 = - =e ml = f ll = f 21 = - = f ml d 2 =e 12 =e 22 = - =e m2 = f 12 = f 22 = "’ = f m2 

d n lrt~^ mn~f ln~f mrv 

> Option 4: columnwise the same design variables (total In design variables) 

d i =eii=e 2 i=-=e ml d 2 =e 12 =e 22 =..=e m2 ... d=e ln =e 2n =.=e mn , ; real parts 

d n + l =fll=f21=-=fml d n+2 = f 12 = f 22 = ‘ ’ ’ = f m2 - d 2n =fln=f2n="=fmn> imaginary parts 


> Option 5: No design variable linking; total 2 mn design variables. 
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Unsteady Aerodynamic Model Tuning 
using Object-Oriented Optimization Tool 


□ 


❖ 


❖ 


The NASA Dryden has developed an Object-Oriented Optimization (O 3 ) tool. 

The O 3 tool leverages existing tools and practices, and allows the easy integration and adoption of 
new state-of-the-art software. 

Local gradient based optimizer as well as 
global optimizers are available. Hybrid methods 
are also available. 

> Optimizers: 

DOT (local), Genetic Algorithm (GA), & 

Big Bang-Big Crunch (BBBC) algorithi 


> Hybrid optimizers: 
GA(CDV)+DOT(CDV), 
GA(CDV)+DOT(CDV)+GA(DDV), 
BBBC(CDV)+DOT(CDV), & 
BBBC(CDV)+DOT(CDV)+BBBC(DDV) 



Frequencies 
Measured 
from Flight 
Test 
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Pre-Processor 
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Flutter Analysis 
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Post-Processor 
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Altitude (ft) Mach 


Time-Invariant Flight Conditions 
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9800 

9700 

9600 
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9200 

9100 

9000 


Time (sec) 
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a> 

■a 

3 


Time (sec) 


10000 

9900 

9800 

9700 

9600 

9500 

9400 

9300 

9200 

9100 

9000 


i 




During Flutter 


Near Flutter 


<■ 


83720 83730 83740 83750 


84660 84670 84680 84690 


Time (sec) 


□ 

First selection: 

Mach = 0.390 

□ 

Second selection: 

Mach = 0.456 

□ Flutter condition: 

Structural Dynamics Group 

Mach = 0.502 


Time (sec) 

Altitude = 9934 ft 
Altitude = 9858 ft 

Altitude = 9837 ft (time varying) 
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Results 



Numerical and measured frequencies (Hz) of the ATW2 during flight test 


Mode 

Natural Frequencies (Hz) 

Measured Aeroelastic Frequencies (Hz) 

Test Validated 
FEM 

GVT 

After take off 
M = 0.390 

Near flutter 
M = 0.456 

During flutter 
M = 0.502 

1 

17.45 

17.45(0.623%) 




2 

43.48 

43.72(0.610%) 

40.45 

38.99 

37.69 

3 

82.98 

83.66(0.778%) 




4 

133.6 

N/A 




5 

153.8 

142.3(0.674%) 





The second aeroelastic frequency before and after unsteady aerodynamic model tuning 


and corresponding scaling factors 


Mach 

Number 

Measured 

(Hz) 

Altitude 

(ft) 

Before Tuning 
(Hz) 

Scaling Factor 
(design variable) 

After Tuning 
(Hz) 

0.390 

40.45 

9934 

41.12 

1.2579 

40.45 

0.456 

38.99 

9858 

40.10 

1.2719 

38.99 
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Altitude (ft) 


Flutter boundaries at Mach=0.502 before and after 
unsteady aerodynamic model tuning 



32000 

30000 

28000 

26000 

24000 

22000 

20000 

18000 

16000 

14000 

12000 

10000 

8000 

6000 

4000 

2000 

0 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


Mach 
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Measured and computed flutter boundaries at Mach = 0.502 


Comment 

Scaling Factor 

(design variable) 

Flutter Speed 

Altitude 

ft 

Flutter Frequency 

Keas 

% difference 

Hz 

% difference 

Measured 

N/A 

276.4 

0.00 

9836.9 

37.69 

0.00 

Before tuning 

1.0 

311.3 

13.0 

3561.5 

37.67 

-0.05 

Use M=0.390 Aero 

1.2579 

277.3 

0.33 

9670.0 

37.69 

0.00 

Use M=0.456 Aero 

1.2719 

276.0 

-0.14 

9912.5 

37.68 

-0.03 
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Conclusions 



□ Model tuning based on flight test data is needed to minimize uncertainties in the unsteady aerodynamic model 
and to increase the safety of flight. 


□ After model tuning (for ATW 2 case) 

❖ Maximum of 13%, 0.13=(311.3-276.4)/276.4, flutter speed error becomes -0.14 %. 


Flutter Boundaries 

Flutter Mach Number 

Altitude = 9836.9 ft 

Flutter Speed 
Differences 

Measured/1.15 = V d 

0.437 

0% 

Measured = 1.15 V d 

0.502 

15% 

Test validated FEM & unsteady aerodynamics; use M=0.456 aerodynamics 

0.5015 

14.8% 

Test validated FEM & unsteady aerodynamics; use M=0.390 aerodynamics 

0.5039 

15.4% 

Test validated FEM; using #2 GVT data; with measured damping 

0.576 

32% 

Test validated FEM; using #2 GVT data; with 3% structural damping 

0.616 

41% 

Test validated FEM; using #1 GVT data; with 3% structural damping 

0.665 

52% 

FEM; before model tuning; with 3% structural damping 

0.762 

74% 


Validated Structural 
Dynamic Model 


Validated Unsteady 
Aerodynamic Model 


Recommended 
Flutter Margins 


r — — - 


Yes 

Yes 

No 


Yes 

No 

No 


15% 

49% 

54% 


ATW2 Case 


14.8 - 15.4 % 
32 - 52 % 
74% 


Structural Dynamic 


Chan-gi Pak-35 


Questions? 



Summary of the Modal Participation Factors 


Mode 

Frequency 

Modal Participation Factor 

Mach 

= 0.60 

Mach 

= 0.75 

Mach 

= 0.82 

Mach 

= 0.95 

1 

17.60 Hz 

68.1 % 

95.5% 

72.9% 

96.2 % 

75.0% 

96.6% 

79.7% 

97.6% 

2 

23.26 Hz 

22.2 % 

18.3% 

16.8 % 

13.6% 

3 

93.99 Hz 

5.2% 

5.0% 

4.8% 

4.3% 

4 

135.4 Hz 

0.0 % 


0.0 % 


0.0 % 


0.0 % 


5 

163.1 Hz 

3.3 % 


2.9 % 


2.6 % 


1.9 % 


6 

174.5 Hz 

0.0 % 


0.0 % 


0.0 % 


0.0 % 


7 

257.5 Hz 

0.7% 

4.5 % 

0.6 % 

3.8 % 

0.5% 

3.4% 

0.3 % 

2.4 % 

8 

391.6 Hz 

0.0 % 


0.0 % 


0.0 % 


0.0 % 


9 

394.3 Hz 

0.1% 


0.1 % 


0.1% 


0.0 % 


10 

445.6 Hz 

0.4% 


0.3 % 


0.3% 


0.2 % 



□ Participation of the first three modes is a function of Mach number. 

□ In-plane modes do not participate for the first flutter mechanism at all. 

❖ Modes 4, 6, and 8 

□ Primary Modes: Modes 1, 2, and 3 

❖ Frequency error should be less than 3%. 

□ Secondary Modes: Modes 4 through 10 (higher) 

❖ Frequency error should be less thanl0%. 
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Flutter Results Before & After Model Tuning 


Mode 

Frequency 

Modal Participation Factors before Model Tuning 

Mach = 0.60 

Mach = 0.75 

Mach = 0.82 

Mach = 0.95 

1 

17.60 Hz 

68.1 % 

95.5 % 

72.9 % 

96.2 % 

75.0 % 

96.6 % 

79.7 % 

97.6 % 

2 

23.26 Hz 

22.2 % 

18.3 % 

16.8 % 

13.6 % 

3 

93.99 Hz 

5.2% 

5.0 % 

4.8% 

4.3 % 

Mode 

Frequency 

Modal Participation Factors after Model Tuning 

Mach = 0.60 

Mach = 0.75 

Mach = 0.82 

Mach = 0.95 

1 

17.45 Hz 

5.0 % 

99.7 % 

8.2% 

99.6 % 

10.7 % 

99.5 % 

22.6 % 

96.4 % 

2 

43.72 Hz 

93.8% 

90.0% 

87.2% 

71.4% 

3 

83.66 Hz 

0.9 % 

1.4% 

1.6 % 

2.4% 


— 

Mach = 0.60 

Mach = 0.75 

Mach = 0.82 

Mach = 0.95 

Before 

Tuning 

Speed 

453.0 Keas 

421.5 Keas 

407.4 Keas 

377.9 Keas 

Frequency 

23.18 Hz 

22.97 Hz 

22.86 Hz 

22.53 Hz 

Altitude 

-7501 ft 

8751 ft 

15010 ft 

25590 ft 

After 

Tuning 

Speed 

337.9 Keas 

340.5 Keas 

341.5 Keas 

344.7 Keas 

Frequency 

35.96 Hz 

35.11 Hz 

34.59 Hz 

32.91 Hz 

Altitude 

8642 ft 

19400 ft 

23475 ft 

29700 ft 

Speed Difference 

33.8 % 

23.8 % 

19.3% 

9.6% 
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